A modified dual-slope method for heat capacity measurements of condensable gases 
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A high resolution non-adiabatic method for measuring the heat capacity (C„) of bulk samples of 
condensable gases in the range 7.5-70 K is described. In this method Cv is evaluated by directly 
comparing the heating and cooling rates of the sample temperature for two algebraically independent 
heat pulse sequences without explicit use of the thermal conductance between sample and thermal 
bath. A fully automated calorimeter for rapid measurement of Cv of molecular solids utilizing this 
technique is presented. 



C/2 



o 



> 

00 

o 

cn 

o 
o 
o 



I 

O 

o 



I. INTRODUCTION 

The adiabatic heat-pulse calorimetry has been widely 
used to investigate the thermodynamic properties of ma- 
terials for more than a century. The applied principle, 
C — AQ/AT, which describes the heat capacity C of 
the sample, is determined by the pulse heat AQ sup- 
plied to the sample under adiabatiG-||Gpnditions and the 
temperature rise AT, is well known.Era Due to the inher- 
ent simplicity as well as the general applicability indepen- 
dent of the sample thermal conductivity, this method is 
the most favored choice for heat capacity measurements 
of condensable gases which have poor thermal conductiv- 
ity in their low temperature solid phase.Em Although the 
traditional adiabatic calorimetry has high precision and 
can be used to determine the latent heat at strong first- 
order transitions, it is very difficult to achieve the resolu- 
tion needed to characterize the temperature dependence 
of Cy{T) (or Cp{T)) close to the critical temperature Tc 
for a second-order transition. Also, because of the inher- 
ent limitations on achieving the ideal adiabatic conditions 
at low temperatures and the long time required to cover 
a few tens of kelvin temperature range with reasonable 
number oiLclaia points, new user friendly non-adiabatic 
techniqueajilj with excellent sensitivity are needed to 
study the heat capacity of condensable gases. Among 
these new techniques, the most sensitive jinethod is the 
ac method devised by Sullivan and SeidelE While the ac 
method allows one to obtain accurate values of C„ , it suf- 
fers from the fact that it normally must be used for small 
samples at low ac frequencies and over limited temper- 
ature ranges (typically T < 20 K). In order to keep the 
sample in equilibrium with the heater and thermometer 
during one cycle, the time constant for thermal relax- 
ation (ti) of the sample plus sample holder (or calorime- 
ter) to the thermal reservoir must be short compared to 
the period of the driving flux. In addition, the samples's 
internal equilibrium time constant {T2) with the sample 
holder (or calorimeter) must be short compared to n . It 
is this latter constraint which restricts the applicability 
of this method to small samples with sufficiently large 
thermal conductivity such as pure metals below 20 K. 

A slightly less sensitive method involving a fixed heat 
input followed by a temperature decay measurement is 



the relaxation method.Q^lljllj In this method the sample 
is raised to an equilibrium temperature above the ther- 
mal reservoir and then allowed to relax to the reservoir 
temperature without heat input. A recent improvement 
of this technique utilizing advanced numerical methods 
was provided by Hwang et aZ.lla The relaxation method 
requires extensive calibration of the heat losses of the 
sample as a function of temperature difference between 
the reservoir and the final sample temperature and re- 
lies on accurate, smooth temperature calibration of ther- 
mometers. The technique depends heavily on the nu- 
merical techniques used to determine many equilibrium 
heat losses during the heating portion of the relaxation 
cycle to obtain good temperature resolution of the heat 
capacity changes. Though Hwang et alx3 addressed the 
problems arising from large r2, this method still fails to 
provide accurate Cy values above '^20 K for large samples 
with poor thermal conductivity. 

A variation of the relaxation method called dual-slope 
method, was first discussed by Riegel and WeberB In this 
method an extremely weak thermal link to the reservoir 
is used and the temperature of the sample is recorded 
over a 10 h. cycle while heating at constant power for 
one half of the cycle, and then allowing the sample relax 
with zero heat input. The heat loss to the reservoir and 
surroundings can be eliminated from the calculation of 
Cv{T) using this technique provided the sample, sample 
holder, and the thermometer are always in equilibrium 
with each other (the reason for long 10 h. cycle to cover 
3 K range), and the reservoir temperature is held con- 
stant over the 10 h. cycle. A further modification of 
this technique is a hybsid between the ac method and 
the dual-slope methodJ3 which reduces the duration of 
the cycle to about two hours. Although the dual-slope 
method is very elegant and easy to implement, the tech- 
nique has only been implemented for small samples (typ- 
ically less than 0.5 g) with good thermal conductivity 
and at temperatures less than ~20 K. The success of this 
method depends heavily on achieving very good ther- 
mal equilibrium between the sample, sample holder, and 
the thermometer, and for large samples with poor ther- 
mal conductivity (i.e., large T2) this method may also fail 
when we consider only the first order approxiniatipn of 
the heat balance equations of Riegel and Weber J3I 
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In the case of condensable gases, the sample size must 
be at least a few grams to reduce the spurious effect 
of sample condensed in the fill line and the heat leak 
through this tube. Due to the low thermal conductiv- 
ity of these samples, in particular for powdered samples 
which may not wet the calorimeter walls, the ti can be 
very large leading to the failure of most of the above 
techniques. The present article describes a modified dual- 
slope method in which C„ is evaluated by directly com- 
paring the heating and cooling rates of the sample tem- 
perature for two algebraically independent heat pulse se- 
quences without explicit use of the thermal conductance 
between sample and thermal bath. For the specific ge- 
ometry of the calorimeter that we used, which is most 
suitable for the heat capacity measurements of condens- 
able gases in the presence of external electric or magnetic 
fields, higher order heat balance equations are calculated, 
which can easily be adapted for other sample configura- 
tions as well. Because of the explicit consideration of the 
higher order equations, the problem of large T2 is natu- 
rally addressed and during the heating and cooling cycles 
one does not require the sample to be in equilibrium with 
the thermometer, or in other words, the thermometer 
need not necessarily record the actual sample tempera- 
ture. Due to this freedom we can obtain Cv of samples 
with poor thermal conductivity very rapidly. 



II. THEORY 

To model the thermal response of a heat-pulse 
calorimeter (Fig. ||) for heat capacity measurements, 
with the schematic diagram shown in Fig. ^ where the 
calorimeter is heated with the application of heater power 
P(i) to some desired temperature T^ax and then allowed 
to cool-down, we start with the following set of heat bal- 
ance equations: 



Pit) = c'fi + A,(r,; - T,) + A,(r^ - Tq 



and 



for heating, and 

= c't'^ + A,(r,' - r,) \r{T[ - To) 

and 

= CTc + A,(T,-T,')-^o(Tc) 



(1) 
(2) 
(3) 
(4) 



for cooling, where C and T are the heat capacity and 
the temperature of the thermometer well plus the inner 
conductor respectively, C and T' are those of the outer 
conductor plus the heater coil. From these equations 
we need to obtain the quantity (C -I- C") which is the 
net heat capacity of the calorimeter. Once the sample 
is condensed in the calorimeter, both C and C will be 



modified but the above equations are still valid. Po{T) 
is the parasitic stray heat due to the conduction of heat 
through the fill line, the stainless steel suspension rod, 
and the thermometer leads as well as radiation from the 
pumping line. From the geometry of the calorimeter (Fig. 
0) it is reasonable to assume that this stray heat affects 
only the inner conductor heat balance equations (Eqs. 2 
& 4). 

Since the thermometer is firmly coupled to the in- 
ner conductor through the thermometer well, the tem- 
perature T is what the thermometer practically records. 
Hence we need to eliminate T' to obtain (C + C). For 
Th = Tc Eqs. 1 & 3 reduce to 

Pit) = C"(T^ - t) + (As + XM - (5) 
From Eqs. 2 & 4 and their first derivatives we obtain 



iU-Z) 



a: 



in - Tc), 



c 



iT,[-T^)^in^T,) + -in-n). 



(6) 



(7) 



After eliminating T^ , and their time derivatives from 
Eqs. 5, 6, and 7 we obtain 

Pit) = iC + C')ifu - t) + ^in - Te) + ^in ~ Te). 

As As 

(8) 
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FIG. 1. Simplified diagram of the thermal system of sam- 
ple, sample holder, and thermal bath. The thermal link be- 
tween the outer conductor plus heater of heat capacity C' at 
temperature T' and temperature-controlled reservoir bath at 
To is represented by the thermal conductance A,.. Ar may be 
controlled by changing the He exchange gas pressure. The 
heater input power P(t) is obtained from a waveform genera- 
tor and a suitable power amplifier. The thermal link (through 
the sample) between the inner conductor plus thermometer 
well of heat capacity C at temperature T and the outer con- 
ductor is represented by the thermal conductance As. 
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Clearly, Eq. (8)|£educes to the first order derivations of 
Riegel and Webem for As 3> 1. However, for condensed 
gases, in particular for powdered samples, this condition 
is never met. In the present geometry, for a sufficiently 
weak thermal link where As ^ A^ , one can ignore the 
last term in Eq. (8). To solve the remaining equation 
exactly for C + C", we need to obtain Eq. (8) for two 
algebraically independent pulse sequences of P{t). 



III. EXPERIMENT 

Figure ^ shows the design of the calorimeter optimized 
for the above technique. The thermal reservoir is made 
up of a brass vacuum can (18 cm long, 4 cm diameter) 
with a manganin wire wound uniformly on the entire 
length of the outer surface. A dilute solution of GE- 
varnish is used to soak the cloth insulation of the man- 
ganin heater wire to provide good thermal contact with 
the vacuum can. A grooved copper tube (4 cm long) is 
soldered to the top flange, which can be used for heat 
sinking the heater and thermometer wires before con- 
necting them to the calorimeter. A 122 cm long stainless 
steel tube (0.96 cm diameter) is soldered to the top flange 
to pump the vacuum space as well as support the vac- 
uum can. copper radiation baffles are placed inside this 
tube at regular intervals to reduce the thermal radiation. 
In addition the top flange supports up to ten individual 
vacuum feed-throughs (not shown), which are thermally 
anchored to the flange through STYCAST 2850FT vac- 
uum seal which has good thermal conductivity at low 
temperatures. The calorimeter consists of two concen- 
tric, thin walled, gold plated, OFHC copper cylinders 
10.2 cm in length providing a 1 mm cylindrical shell space 
for condensing the samples. This shell space is vacuum 
sealed at both ends with, homemade 1 mm thick STY- 
CAST 2850FT O-ringsJia The outer conductor (2.7 cm 
inner diameter and 0.2 mm wall thickness) with a brass 
heater wire uniformly wound on the entire length of its 
outer surface serves as a radial heating element. Uniform 
winding of the heater wire is very crucial in eliminating 
longitudinal heat flow thereby reducing the temperature 
gradient along the length of the calorimeter. For good 
thermal contact with the outer cylinder, a dilute solu- 
tion of GE-varnish was applied to the heater wire. In 
the middle of the inner conductor (2.5 cm diameter, 0.4 
mm wall thickness) a thin copper ring as well as a copper 
thermometer well are soldered (see Fig. 0). The larger 
wall thickness (0.4 mm instead of 0.2 mrn) for the inner 
cylinder is necessary to withstand the pressures gener- 
ated when the solid samples melt at high temperatures. 
The copper ring sandwiched between two threaded nylon 
wings supports the calorimeter when suspended from the 
top flange with the help of a threaded stainless steel rod. 
This arrangement not only reduces the unwanted heat 
leaks to the calorimeter but also provides excellent me- 
chanical stability when the free end of the stainless steel 



rod is slipped into the post at the bottom of the vacuum 
can (see Fig. ||). 
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FIG. 2. Cross section of the thermal reservoir and calorime- 
ter, (a) Top flange with vacuum feed-throughs, (b) brass vac- 
uum can with highly uniform heater wire wound on the outer 
surface, (c) 1.6 cm diameter copper tube with grooves on the 
outer surface, (d) STYCAST 2850FT sealji3 (e) outer con- 
ductor with highly uniform heater wire wound on the outer 
surface, (f)inner conductor, (g) copper thermometer wells sol- 
dered to the walls, (h) threaded stainless steel rod, (i) 0.5 mm 
thick copper ring soldered to the inner conductor, (j) threaded 
nylon wings, (k) sample fill line, (1) copper radiation baffles, 
(m) stainless steel tube for support as well as pumping. 

The 1.6 mm diameter cupronickel alloy tube serving as 
the sample fill line is coiled into a three turn spring (not 
shown) before soldering it to the inner cylinder. To be 
able to apply electric field across the sample cell when 
necessary, it is important to ground the outer conductor 
(to reduce electrical interference from the heater current) 
and apply the high voltage to the inner conductor. Hence 
the coiled section of the fill line is isolated electrically as 
well as thermally from the remaining length with the help 
of a vacuum joint made up of a wider cupronickel tube, 
cigarette paper, and STYCAST seal. To assist sample 
condensation, the manganin heater wire is wrapped on 
the fill line tube outside of the vacuum can. Two, factory 
calibrated germanium resistance thermometeralJ placed 
inside the thermometer wells are used for recording the 
temperatures of the calorimeter and the thermal reser- 
voir. The entire assembly is placed in a second vacuum 
can (125 cm long, 5 cm diameter) which can be lowered 
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into a commercial liquid He dewar. With this arrange- 
ment including a few suitably placed copper radiation 
baffles at the neck of the outer vacuum can, a standard 
100 liter dewar lasted for almost a month while heat ca- 
pacity experiments were carried out in the 4.2 — 75 K 
range. 

IV. RESULTS AND DISCUSSION 

Our observations indicate that due to the high degree 
of cylindrical symmetry, even though the calorimeter is 
rather unusually long (~10 cm), thermometers record 
the expected temperature values to very high accuracy. 
Figure |^ shows typical temperature excursions when the 
calorimeter is filled with pure N2 as a sample. For sim- 
plicity we chose square and triangular voltage pulses 
for two algebraically independent P{t) sequences. Both 
pulses can be easily generated with the help of a standard 
signal generator or data acquisition card and a suitable 
operational amplifier. Before the heat pulse is applied, 
the calorimeter is left to equilibrate with the surround- 
ings until its temperature trace is horizontal with the 
time axis. Typically, when the reservoir temperature is 
adjusted to a new value, the calorimeter attained equi- 
librium in 30 to 60 min. To ensure that the calorimeter 
always attained equilibrium from either a higher temper- 
ature or a lower temperature for all pulse sequences (i.e., 
T — > Tq*" or T — > Tg"), and that small hysteresis of the 
the thermometer does not affect the trace, when Tq is set 
to a higher value, a small heat pulse is simultaneously 
applied to the calorimeter so that its temperature rises 
above the equilibrium temperature by at least 2 to 3 K. 
When such a precaution was not taken, we observed that 
the base line of the two traces in Fig. ^ differed by ~100 
mK resulting in spurious results when fit to Eq. (8). 

Although we assumed that the last term in Eq. (8) 
will be neghgible for As » Ar, it is clear that for the long 
exponential decay section of the traces (below the lower 
horizontal line in Fig. ^) where {Th—Tc) can be very large, 
the above approximation may not be valid. We observed 
that for a typical 6.5 K pulse (similar to Fig. ^, the useful 
window of temperature where the above approximation 
is valid for obtaining faithful results is only 3.5^4 K (the 
section between the two dashed lines in Fig. ^). It is cru- 
cial that during the entire pulse sequence, the reservoir 
temperature To should be conjStant. With the help of a 
program written in Lab VIEW Jla which controls the volt- 
age across the-, reservoir heater through a multifunction 
AD/DA carets and with the thermometer output in the 
negative feedback loop, we were able to maintain Tq to 
within 30 mK during the entire pulse sequence (see Fig. 
^). Because the feedback is controlled through software, 
excellent stability in Tq was achieved by changing the 
feedback parameters depending on the pulse height and 
temperature range. 

Figure |4| shows the heat capacity curves of the 



calorimeter filled with pure N2 for two Tq values. Ex- 
cellent overlap of the curves demonstrates that with this 
technique, we can obtain heat capacity of samples with 
poor thermal conductivity very rapidly (~3.5 K range 
in 3 h.), without further curve fitting generally required 
with other techniques.EI By extending the Lab VIEW pro- 
gram written for controlling the reservoir temperature 
To, we fully automated the data acquisition process to 
obtain pulse sequence data for various Tq values with an 
interval of 2.5 K to ensure data overlap from adjacent 
pulse sequences. Once the data are obtained, standard 
third degree polynomial fits (only for the useful temper- 
ature interval similar to the one between dashed lines in 
Fig. H) are employed for obtaining time derivatives in Eq. 
(8). Once the heat capacity is calculated for each pulse 
sequence, no further curve fitting is necessary and in gen- 
eral for r < 30 K, the overlap between adjacent curves 
is better than the one shown in Fig. ^. Figure || shows 
the heat capacity of the empty calorimeter (with small 
amounts of He gas for better thermal contact) obtained 
through the above automated process in the 10 — 70 K 
range. The kink at 60 K, the presence of which is acciden- 
tal, shows the heat capacity near a continuous phase tran- 
sition with excellent resolution. The data shown without 
further curve fitting clearly demonstrates the high sensi- 
tivity of the technique. 
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FIG. 3. Traces of T{t) and To(t) for two heat pulses when 
the calorimeter is filled with pure N2 sample. Useful interval 
for obtaining heat capacity is shown between two dashed lines. 
The difi'erence between the trace of To and the baseline of T 
is due to To- 
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In summary, we have demonstrated the use of a novel 
technique to study the heat capacity of moderately large 
samples with poor thermal conductivity in the 7.5-70 K 
range. A fully automated calorimeter for rapid measure- 
ment of the heat capacity of condensable gases utilizing 
the above technique has been presented. The technique 
along with the automated calorimeter with a provision 
to apply external electric and magnetic fields is particu- 
larly useful for the study of continuous phase transitions 
in molecular solids as well as field induced changes in the 
heat capacity. 
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FIG. 4. Heat capacity of the calorimeter filled with pure 
N2 sample. Data are obtained for two pulse sequences with 
respective To values separated by 2 K. 
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FIG. 5. Heat capacity of the empty cell with small amounts 
of He exchange gas for better thermal contact. No curve fit- 
ting techniques are utilized as the adjacent heat pulse data 
overlap very well. The sharp kink at 60 K, whose origin is 
not known, clearly demonstrates the high resolution of the 
technique near continuous phase transitions. 
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